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The A decay rates in nuclei has been calculated in ref. El using the Propagator 
Method in Local Density Approximation. We have studied the dependence of the 
widths (including the one for the two-body induced process ANN — > NNN) on 
the NN and AA short range correlations. Using a reasonable parametrization of 
these correlations, as well as realistic nuclear densities and A wave functions, we 
reproduce, for the first time, the experimental non-mesonic widths from medium 
to heavy hypernuclei. 

1 Introduction 

The study of hypernuclear physics may help in understanding some present pro- 
blems related, for instance, to aspects of weak interactions in nuclei, or to the 
origin of the spin-orbit interaction in nuclei. Besides, it is a good instrument to 
study the role of quark degrees of freedom in the hadron-hadron interactions at 
short distances and the renormalization properties of pions in the nuclear medium. 

The weak decay of hypernuclei occurs via two channels; the so called mesonic 
channel: 

A — > ttN (T m ), (1) 

and the non-mesonic one, in which the pion emitted from the weak hadronic vertex 
is absorbed by one or more nucleons in the medium, for example 

AN ^ NN (Ti), (2) 

ANN — > NNN (r 2 ). (3) 

Obviously, the non-mesonic processes (IVm = Ii + T 2 + ...) can also be mediated 
by the exchange of more massive mesons. The non-mesonic decay is only possible 
in nuclei and, nowadays, the study of the hypernuclear decay is the only practical 
way to get information on the weak process AN — ► NN, especially on its parity 
conserving part. 

The free A decay is compatible with the AI =1/2 isospin rule, which is also 
valid in other non-leptonic strangeness changing processes. This rule is based on the 
experimental observation that F^^^-p/FA^nOn — 2, but it is not yet understood 
on theoretical grounds. From theoretical calculations like the one in ref. 13 and 
from experimental measurements □ there is some evidence that the AI =1/2 rule 
is broken in the nuclear mesonic decay. However, this is essentially due to shell 
effects and might not be directly related to the weak process. A recent estimate of 
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AI = 3/2 contributions to the A A" — > NN reaction □ found moderate effects on the 
hypernuclear decay rates. In the present calculation of the decay rates in nuclei we 
will assume this rule as valid. The momentum of the final nucleon in A — > -kN is 
about 100 MeV for A at rest, so this process is suppressed by the Pauli principle, 
particularly in heavy nuclei. It is strictly forbidden in infinite nuclear matter (where 
kp ~ 270 MeV), but in finite nuclei it can occur because of three important effects: 
I) in nuclei the Pauli blocking is less effective and the hyperon has a momentum 
distribution which allows larger momenta for the final nucleon, 2) the final pion 
feels an attraction by the medium such that for fixed momentum it has an energy 
smaller than the free one and consequently, the final nucleon again has more chance 
to come out above the Fermi surface, and 3) on the nuclear surface the local Fermi 
momentum is smaller than k F and favours the decay. Nevertheless, the mesonic 
width decreases fast as the mass number A of the hypernucleus increases □. The 
mesonic rate is very sensitive to the pion self-energy, so that from the study of 
the mesonic decays it could be possible to extract important information on the 
pion-nucleus optical potential, which we do not know, nowadays, in a complete 
form. 

The final nucleons in the non-mesonic process AA~ — > AW come out with a 
momentum ~ 420 MeV, so that this decay is not forbidden by the Pauli principle. 
On the contrary, apart from the s-shell hypernuclci, it dominates over the mesonic 
decay. The non-mesonic channel is characterized by large momentum transfers, so 
that the details of the nuclear structure do not have a substantial influence while 
the AW and AN short range correlations turn out to be very important. There 
appears to be an anticorrelation between mesonic and non-mesonic decay modes: 
indeed the total lifetime is fairly constant from light to heavy hypernuclei &u: r exp = 
(0.5 rjr/vee- 

Nowadays, the main problem concerning the weak decay rates is to reproduce 
the experimental values for the ratio r„/r p between the neutron and the proton 



induced widths An — > nn and Ap — > np. All theoretical ca 
the experimental data for all the considered hypernuclei 13' 
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In the One Pion Exchange (OPE) approximation the values for this ratio are 0.1 -r 
0.2. On the other hand, the OPE model has been able to reproduce the 1- 
stimulated non-mesonic rates Ti = T n + T p for light and medium hypernuclei I 
In order to solve this puzzle many attempts have been made up to now, but without 
success. Among these we recall tJ^inclusion in the AA~ — > AW transition potential 
of mesons heavier than the pion EJtffl, the inclusion of contributions that violate the 
AI =1/2 rule u and the description of the short range baryon-baryon interaction 
in terms of quark degrees of freedom 0. This last calculation is the only one which 
has found a fairly large (but not sufficient) increase of the neutron to proton ratio 
with respect to the OPE one. However, this calculation is carried out only for 
s-shell hypernuclei; moreover, the employed quark-lagrangian does not reproduce 
the experimental ratio between the AI = 1/2 and AI = 3/2 transition amplitudes 
for the A free decay. 
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The analysis of the ratio T n /T p is influenced by the 2-nucleon induced process 
ANN — ► NNN. By assuming that the meson produced in the weak vertex is mainly 
absorbed by a strongly correlated neutron-proton pair, the 3-body process turns out 
to be Anp — > nnp, so that a considerable fraction of the measured neutrons could 
come from this channel and not only from the An — ► nn and Ap — > np ones. In 
this way it might be possible to explain the large experimental T n /T p ratio, which 
originally has been analyzed without taking into account the 2-body. .stimulated 
process. Nevertheless, the situation is far from being clear and simple a't3. In fact, 
Tn/r^ is sensitive to the energy spectra of the emitted nucleons, whose calculation 
also_Eequires a careful treatment of the final state interactions of the nucleons. In 
ref. IB the energy distributions were calculated using a Monte Carlo simulation to 
describe the final state interactions. A direct comparison of those spectra with 
the experimental ones seems to favour large r n /T p ratios, in strong disagreement 
with the OPE predictions. However, it was also pointed out the convenience of 
measuring the number of protons per decay event. In fact, this observable gives 
a more reliable r„/r p ratio and it is less sensitive to details of the Monte Carlo 
simulation accounting for the final shape of the spectra. 

2 Decay widths 

We present here the results of a new evaluation of the decay rates for medium to 
heavy A-hypernuclei based on the Propagator Method introduced in ref. 113. The 
parameters of the model were adjusted to reproduce the non-mesonic width of A 2 C. 
Then the decay rates for heavier hypernuclei have been predicted and compared, 
when possible, with the experiment. 

2.1 The model 

The propagator technique, estensively explained in refs. 00, provides a unified 
picture of the different decay channels and it is equivalent to the standard Wave 
Function Method 113. The calculation of the widths is performed in nuclear matter, 
and then extended to finite nuclei via the Local Density Approximation (LDA) . 

The starting point for the calculation is the A — > irN effective lagrangian, 
parametrized in the form: 

Ca-kN = Gm^i\) N (A + Bj 5 )t ■ (fr^A + h.c, (5) 

where the values of the weak coupling constants G, A, B are fixed on the observables 
of the free A decay. In order to enforce the AI =1/2 rule, in eq. (||) the hyperon 
is assumed to be an isospin spurion with I z = —1/2. The decay rate of a A in 
nuclear matter is obtained from the imaginary part of the A self-energy through 
the relation: 

T A = -2ImS A . (6) 

The A self-energy, Ea, contains an infinite series of diagrams in which the pion 
emitted in the A vertex couples to lp-lh, 2p-2h, etc excitations, giving rise to the 
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different decay channeis when they are cut in all the possible ways. The structure 
of the decay width in nuclear matter is the following: 



G ^l) - 9 x2 ~ toTI ' ( 8 ) 



r A (fc) oc J dq...6{\ k-q\ -k F )e(ko-E N (k-q)-V N )ImG„(q) \ qo=ko - Elf tf-<r)-v N > 

(7) 

where 

1 

Qo ~ f ~ ml - T,*(q) : 

is the pion propagators in nuclear matter, E* being its proper self-energy; hp is 
the Fermi momentum, En is the nucleon total free energy and Vjv is the nucleon 
binding energy. Moreover, k = E\(k) + Va, is the A energy, containing a binding 
term (which we take from the experiment). 

The decay widths in finite nuclei are obtained in LDA; a local Fermi sea of 
nucleons, related to the nuclear density by 



3 



1/3 



fcF(f) = ||7rV(r)| , (9) 

is introduced. Besides, the nucleon binding potential Vn also becomes r-dependent 
in LDA. The decay width in finite nuclei is then obtained through the relations: 



r A (fc) = I df\ ^ A (r) | 2 r A [fc, P (r)}, r A = J dk\ </> A (fc) | 2 r A (fc), (10) 

where if)\ ("0a) is the A wave function in the coordinate (momentum) space. 



2.2 Results 

In ref. we have studied the influence of the short range correlations and the A 
wave function on the decay width of A 2 C, which we have used as a testing ground 
for the theoretical framework in order to fix the parameters of our model. 

The baryon-baryon correlations are governed by the Landau parameters g' and 
g' A . The choice of these parameters which reproduce the experimental non-mesonic 
rate of A 2 C is g' = 0.8, g' A = 0.4. Neglecting the contribution of the 2-body induced 
decay channel (namely in ring approximation), with the same g A , we find that 
g' = 0.7 should de used. This-walue is in agreement with the fenomenology of the 
(e, e') quasi-elastic scattering til. 

The A wave function is obtained from a Wood-Saxon well which exactly repro- 
duces the first two experimental A single particle eigenvalues (s and p levels). The 
numerical results for \ 2 C are^hown in table.f], .where they are compared with the 
experimental data from BNL ll3 and KEK In addition to the Wood-Saxon 

wave function that reproduces the s and p levels, we also made a test using the 
harmonic oscillator wave function with an "experimental" frequency uj obtained 
from the s-p experimental energy shift, the Wood-Saxon wave function of Dover 
et al. Ej and the microscopic wave funclipn calculated from a non-local self-energy 
using a realistic AN interaction in ref. E3. We have found that all these different A 
wave functions give rise to total decay widths which may differ at most by 15%. 
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Table 1. Decay rates for A 2 C 





Calculation 


BNL ta 


KEKB 


KEK New BE3 


rv 


0.25 


0.11 ±0.27 


0.36 ±0.15 


> 0.11 




0.82 








r 2 


0.16 








TjVM 


0.98 


1.14 ±0.20 


0.89 ±0.18 


< 1.03 


r T 


1.23 


1.25 ±0.18 


1.25 ±0.18 


1.14 ±0.08 



Table 2. Decay rates of medium-heavy A-hypernuclei 



A Z 




Ti 


r 2 


r T 


12 p 
A ^ 


0.25 


0.82 


0.16 


1.23 


|«Si 


0.07 


1.02 


0.21 


1.30 


l°Ca 


0.03 


1.05 


0.21 


1.29 


fFe 


0.01 


1.12 


0.21 


1.35 


89v 
A Y 


6 • 10~ 3 


1.16 


0.22 


1.38 


A 39 La 


6 • lO- 3 


1.14 


0.18 


1.33 


A 08 Pb 


1 • 10" 4 


1.21 


0.19 


1.40 



Then, using Wood-Saxon wave functions that reproduce the s and p A-levels 
and the Landau parameters g' = 0.8, <? A = 0.4, we have extended the calculation to 
heavier hypernuclei. We note that, in order to reproduce the experimental s and p 
levels for the hyperon we must use potentials with a nearly constant depth, around 
28 -7- 32 MeV, from medium to heavy hypernuclei. 

Our results are shown in table ^. The mesonic rate rapidly vanishes by increasing 
the mass number A. This is well known and it is related to the decreasing phase 
space allowed for the mesonic channel, and to smaller overlaps between the A wave 
function ip\ and the nuclear surface, as A increases. The 2-body induced decay is 
rather independent of the hypernuclear dimension and it is about 15% of the total 
width. The total width is also nearly constant with A, as we already know from the 
experiment. In fig. |l| we compare the result^cara-tahle || with recent (after 1990) 
experimental data for non-mesonic decay ll3li3 ! t3Eill23 ! E2l. We remind that the data 
for nuclei heavier than A 2 C refer to the total width. The theoretical results are in 
good agreement with the data over the whole hypernuclear mass range explored. 
Moreover, we also see that the saturation of the AN — * NN interaction in nuclei 
appears to be well reproduced. 

The main open problem in the study of weak hypernuclear decays is to under- 
stand the large experimental value of the ratio T n /T p . However, we have to remind 
that the data for T n /T p have a large uncertainty and they have been analyzed with- 
out taking into account the 3-body decay mechanism. The study of ref. E2l showed 
that, even if the three body reaction is only about 15 % of the total decay rate, this 
mechanism influences the analysis of the data determining the ratio r„/r p . Energy 
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spectra of neutrons and protons from the non-mesonic decay mechanisms were cal- 
culated in ref. E_j. The momentum distributions of the primary nucleons were 
determined from the Propagator Method and a subsequent Monte Carlo simulation 
was used to account for the final state interactions. It was shown that the shape 
of the proton spectrum is sensitive to the ratio T n /T p . In fact, the protons from 
the three-nucleon mechanism appeared mainly at low energies, while those from 
the two-body process peaked around 75 MeV. Since the experimental spectra show 
a fair amount of protons in the low energy region they would favour a relatively 
large three-body decay rate or, conversely, a reduced number of protons from the 
two-body process. Consequently, the experimental spectra were compatible with 
values for T n /T p around 2 3, in strong contradiction with the present theories. 

Within the Propagator Method with modified parameters and realistic A wave 
functions we have then generated nest nucleon spectra for various hypernuclei using 
the Monte Carlo simulation of ref. c3. Although the new non-mesonic widths are 
sizably reduced (by about 35%) with respect to those of refs. IMM, the resulting 
nucleon spectra are practically identical. The reason is that the ratio r 2 /Pi is es- 
sentially the same in both models, and the momentum distributions for the primary 
emitted nucleons are also very similar. As a consequence, the conclusions drawn in 
ref. B still hold. 

Therefore, the origin of the discrepancy between theory and experiment for 
the ratio r n /r p still needs to be resolved. However, we must notice that on the 
experimental side, although new spectra are now available LLj'EZl, they have not been 
corrected for energy losses inside the target and detector, so a direct comparison 
with the theoretical predictions is not yet possible. Attempts to incorporate these 
corrections by combining a theoretical model for the nucleon rescattering in the 
nucleus with a.-simulation of the energy losses in the experimental set-up are now 
being pursued Ej. On the other hand, a forward step towards a clean extraction of 
the ratio T n /T p would be obtained if the nucleons which come out from the different 
non-mesonic processes, AA" — > AW and ANN — » NNN were disentangled E3. 
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